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Abstract 

We consider the dynamics of timelike spherical thin matter shehs in vacuum. A general formahsm 
for thin shells matching two arbitrary spherical spacetimes is derived, and subsequently specialized 
to the vacuum case. We first examine the relative motion of two dust shells by focusing on the 
dynamics of the exterior shell, whereby the problem is reduced to that of a single shell with 
different active Schwarzschild masses on each side. We then examine the dynamics of shells with 
non- vanishing tangential pressure p, and show that there are no stable — stationary, or otherwise — 
solutions for configurations with a strictly linear barotropic equation of state, p = aa^ where a is 
the proper surface energy density and a G (—1,1). For arbitrary equations of state, we show that, 
provided the weak energy condition holds, the strong energy condition is necessary and sufficient 
for stability. We examine in detail the formation of trapped surfaces, and show explicitly that 
a thin boundary layer causes the apparent horizon to evolve discontinuously. Finally, we derive 
an analytical (necessary and sufficient) condition for neighboring shells to cross, and compare the 
discrete shell model with the well-known continuous Lemaitre-Tolman-Bondi dust case. 
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I. INTRODUCTION 



The matching of two arbitrary spacetimes along a given hypersurface plays an important 
role in general relativity, with a rich plethora of applications, such as the dynamics of thin 
matter shells construction of cosmological models 0, collapse of bounded bodies 0, and 
wormholes 0. The standard techniques to achieve such matching are not circumscribed to 
four-dimensional manifolds, and can be easily applied to higher- dimensional cases, such as 
n-branes in the so-called brane world cosmology [Q. 

Within the context of classical general relativity, the thin shell matching problem was 
first studied by Sen ||^, Lanczos [0] and Darmois [Q, and later further developed by Israel 0, 
who produced a coordinate invariant formalism by applying the Gauss-Codacci equations to 
a non-null three-dimensional hypersurface imbedded in a four-dimensional spacetime. The 



null case was studied by Barrabes and Israel ||l0[- On the shell, Einstein's equations reduce 
to the Lanczos equation, where the jump of the extrinsic curvature across the shell plays 
the role of the four- dimensional Ricci tensor, being thus related to the surface stress-energy 
tensor. Comprehensive reviews of the matching problem in general relativity may be found 

in og. 

In this paper, we first consider two shells in vacuum, and analyze the dynamics from 
the viewpoint of the exterior shell, thereby reducing the problem to that of a single shell 
immersed in two adjacent Schwarzschild spacetimes with different mass parameters on each 
side. For shells with pressure, we show that (i) a strictly linear barotropic equation of state 
is incompatible with stable (oscillatory or stationary) solutions, and (ii) for matter with a 



real local sound speed {Cg = ^/dp/da > 0) obeying the weak energy condition, the strong 
energy condition is necessary for stability, regardless of the equation of state. We then study 
the formation of trapped surfaces in the spacetime, and show explicitly that the existence 
of a thin matter shell introduces a discontinuity in the apparent horizon curve. The related 
problem of shell crossings is examined, and an analytical necessary and sufficient condition 
for neighboring shells (defined in a precise manner) to cross is derived. We further show 
that, although shell crossing singularities occur in a multi-dust-shell case just as they do in 
the continuous Lemaitre-Tolman-Bondi (LTB) dust case, the former cannot be taken as the 
discrete analogue of the latter, since, e.g., the energy density remains finite at discrete shell 
crossings, whereas it diverges (together with some curvature components) in the LTB case. 
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Physically, this is related to the fact that individual shells move geodesically in the LTB 
spacetime (i.e., proper time equals comoving time, for each shell), but the same does fails 
to hold in the discrete case: dust particles do move geodesically on each shell, but, from the 
viewpoint of the four-dimensional spacetime, they are accelerated. 

The paper is organized as follows: Section II derives a general formalism for the matching 
of two arbitrary spherical spacetimes across a timelike thin shell with arbitrary matter 
content. In Sec. Ill we specialize to the case of two dust shells in vacuum; the general 
relativistic equation of motion for the exterior shell is derived, as well as its Newtonian 
analogue. The well-known case of a single dust shell in vacuum is trivially recovered. In 
Sec. IV, the dynamics of a dust shell with a LTB interior is discussed, and the boundary layer 
case of closed Lemaitrc-Fricdmann-Robertson- Walker (LFRW) models is readily obtained as 
special case. Section V studies shells with a non-vanishing tangential pressure component, 
and several stability results are produced. In Sec. VI, we derive the condition for formation of 
apparent horizons in an arbitrary spherical spacetime, and then examine the particular case 
of a dust shell collapsing in Schwarzschild background. Section VII discusses the occurrence 
of shell crossing in a spacetime with a finite number of shells, derives an explicit condition 
for shell crossing of neighboring shells, and compares it with the continuous dust LTB case. 
Section VIII concludes with a summary and discussion. 

Natural geometrized units, in which G — c— 1, are used throughout. 



II. MATCHING FORMALISM FOR ARBITRARY SPHERICAL SPACETIMES 

The complete four-dimensional spacetime consists of an interior region connected 
to an exterior Ai^ by a timelike three-dimensional thin shell S. The Ai± regions are 
characterized by the spherical metric: 

dsl = -hldtl + fldr^ + r'(ifi^ (2.1) 

where /i±, f± are functions of t± and r, and the coordinate systems are adopted. On 

E there is a natural holonomic basis {e(a)} given by 

= (2-2) 
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where j^*^, a = 0, 1, 2} are intrinsic coordinates on S. The induced three-metric 'jab on S is 
then 

lab = 9t,ue'"^af'(b), (2-3) 

and it is the same on both sides of S, since the four- metric must be continuous across it. 
The surface S is parametrically defined by 

= r-rE(t±) =0. (2.4) 

Continuity of g^^, across S requires rs(t±) = -R(t), whence 

dsl = -dr^ + R\T)dQ\ (2.5) 
^ = h-^'^l + flR^^rj^, (2.6) 

where R = dR/dr, and r is the proper time measured by a comoving observer with the 
shell, with four- velocity 

= v±5l + R5^. (2.7) 
The spacelike (outward pointing) unit normal to S is 

< = 9^'a-'^^^ = /i±V±i?C + fi'h^v±Si^, (2.8) 

where = n^rv^, and n^u^ = 0. We shall henceforth drop the ± subscript for clarity and 
without detriment, since all the formulae will apply equally to both sides of S. 
The normal extrinsic curvature, Kab, is |0 

Kab := -^Me(fe)V.ef„) = -n^ [g^b + ^^^-Q^Q^ ) ' ^^.9) 

which is related to the (yet unspecified) surface stress-energy tensor Sab on E via the Lanczos 
equation 0] 

[Kab]-lab[K] = -8nSab ^ [Kab] = StT {Sab - ^S^ab) , (2.10) 

where 

[Kab] ^ K^, ~ K-, [K]^r'[Kab], S^JabS'"'. (2.11) 
The four-dimensional stress-energy tensor associated with S can be written as a distribution 



as 



T- = 5'^^eJ,)e^,)|a|<5($). (2.12) 
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The Lanczos equation (for the shell), together with the (exterior) Einstein equations, lead 
to the standard Israel junction conditions 0: 

S'^^K^,} ^ \{K^, + K-^S-' = [T,^n'^n% (2.13) 
= [ef„)T,.n1. (2.14) 

These two equations are identities that must be satisfied throughout the time development 
of the shell, and not genuine dynamical equations, since they follow from the momentum and 
Hamiltonian constraints imposed on S the dynamics is encoded in the Lanczos equation 

With the metric (|2.1| ), the non- vanishing components of K^b are: 



Kee = sm~''eK^^ = ^, (2.16) 
where D± = n^V^ is the normal derivative in the n'^ direction, and ' = dr- 



ill. TWO COLLAPSING SCHWARZSCHILD SHELLS 

We shall examine here the case of two neighboring thin shells of dust. Specifically, the 
model consists of an inner shell with gravitational mass m_ and an outer shell with m+, in 
an otherwise empty spacetime. An immediate consequence of Birkhoff 's theorem is that the 
two-shell case can be reduced to that of a single shell with appropriate adjoint spacetime 
metrics. If one is interested in the dynamics of the inner shell, the problem is that of a 
single shell with a Minkowski interior (to ensure regularity of the metric at the center) and a 
Schwarzschild exterior with mass M = m_ +m+. Focusing on the exterior shell, the problem 
reduces to that of a single shell with an interior Schwarzschild metric with mass m_ and 
Schwarzschild exterior with mass M. We shall henceforth adopt this latter viewpoint. The 
metrics m. M.± are then: 

.2 _ 7 2 1.2 , 7,-2 7 2 , „2 7r,2 



^4 = -h^dt' + hl'dr' + r'dQ.\ (3.1) 



with 



2m_ 2M , , 

h^ = \l , h+ = \l , M = m^ + m+. (3.2) 



In this case, Eqs. (|2.15|) - (|2.16|) simplify to 



Krr = -^^, Kee = Rh\ (3.3) 



Now, let a be the total rest mass of the shell per unit proper area. The stress-energy tensor 
Sab on S is then 

Sab = 0-UaUb = ae^^^^~^e\f^-^u^Uy = a6l6l, (3.4) 
and the total rest mass of the shell is 

ms = y Sabu'^u'' ^/det^de A d0 = AnR^a. (3.5) 

The conservation equation ( |2.14| ) gives 



- + 2i^ = T^ = 0, 3.6 

i.e., the proper mass of the shell is conserved during the evolution; microscopically, this 
simply reflects the fact that the total number of particles in the shell is conserved, and thus 
the flux three-vector = cru" is divergence-free. 

The dynamical evolution of the shell can be easily obtained from the 99 component of 
the Lanczos equation [a straightforward calculation shows that Eq. ( p.l3|) is equivalent to 



[Kqq] = —AiraR'^, and the rr component of the Lanczos equation does not yield additional 



information — its first integral is automatically satisfied by Eq. (|3.8| )] 



[Kee] = -AiraR'. (3.7) 
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From Eqs. (|3.2|), (|3.5| ), and ( |3.7| ) we obtain 
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R"=[z^] -1 + ^^ + 7^- (3-8) 



M + m_ m% 



R 4i?2- 



This is the full general relativistic equation governing the motion of the shell S with active 
gravitational mass m_|_ in an interior Schwarzschild background with active gravitational 
mass m_. Since both the gravitational and rest masses of the shell remain constant during 
the evolution, we can introduce the dimensionless constant k = m+/mY., and rewrite the 
dynamical equation solely in terms of the gravitational masses: 

^2 _ p _ M + m. 



Based on its Newtonian analogue (derived below), one may interpret the terms in this 
equation as follows: A;^ — 1 is the total specific (i.e., per unit gravitational mass m^) binding 
energy of the system, the second term is the potential energy, and the third term is a self- 
binding energy for S, which is a relativistic correction to the Newtonian case. We note 
that, since energy (potential, or otherwise) gravitates, the proper mass of the shell will only 
coincide with its active gravitational mass if the total binding energy of the system vanishes, 
when k = 1. For given initial data {ri; -R(ri), i?(ri)}, collapsing or expanding solutions may 
be obtained, depending on the value of G (0, +00). For k G (0,1), potential (negative) 
energy dominates, and there is a maximum radius at which the shell is momentarily at 
rest before recollapsing, i.e., the system is gravitationally bound; for k = 1, the kinetic 
energy vanishes exactly at spatial infinity, and the system is said to be marginally bound; 
if A; G (l,+oo), the kinetic energy is positive-definite at infinity [it "equals" m+(/c^ — l)/2] 
and the system is said to be gravitationally unbound. 



Newtonian limit 



The Newtonian limit is obtained by the requirement of non-relativistic velocities, \R\ -C 1, 
which implies [|13| the weak-field limit, e = m±/R << 1. To first-order in e, Eq. (|3.9| ) reads 



R^ = k^-1 + ^+J^- + o{e'). (3.10) 
R 

To check that this is indeed the correct limit, we present below a simple Newtonian derivation 
of the same equation. Let m^p be the rest masses of the inner and outer shells, respectively. 
We want to follow the motion of the m_|_ shell in the gravitational potential of the m_ shell. 
By Gauss's theorem, this reduces to the problem of finding the dynamical equation of a 
spherical shell S with radius r = R and mass m+, subject to the gravitational potential of 
a point mass m_ located at r = 0. The action for such system is 

S = j L{R,R,t)dt = j [T -{U^ + U^)]dt, (3.11) 

where 

T = ^m^R\ (3.12) 

U- = (3.13) 
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= 2n <^^(r)p{rydr = --±, (3.14) 
Jo '^^ 

are the kinetic, potential, and self-binding (of the shell) energies, respectively. The function 

$£(r) = —m^/r is the self-gravitational potential of the shell, and p(r) = o"5(r — R) is the 

energy density, where a = m+/(47ri?^) is the shell's surface energy density. Extremizing 

6S = yields the Euler-Lagrange equation 

fl + ^^^^ = 0. (3,15) 



which integrates to 



2 ^ , M + m. 



R' = n + , (3.16) 



where M = m+ + m_, and Q is an integration constant, which represents (twice) the total 
binding energy of the system per unit shell mass. This agrees thus with Eq. ( p.lO| ) upon 
setting Q = k'^ — 1. 



B. Single dust shell in vacuum 

By setting m_ = 0, we readily obtain the case of a single spherical dust shell in an 
otherwise vacuum spacetime, originally studied by Israel [0]: 

From the first equality above it follows that 2kR{l — k) < M. If A; > 1, this is trivially 
satisfied, but for k < 1 this imposes an upper bound on R: 

As before, if the system is gravitationally bound (A; < 1), any initially expanding shell will 
reach a maximum radius — uniquely determined by the gravitational mass of the shell and 
its total binding energy — and then collapse back to i? = 0. Gravitationally unbound shells 
(A; > 1) will expand to i? cxo. As noted by Israel, the upper bound ( p.l8| ) implies that 
-Rmax > 2M, and hence there are no timelike stationary shells with R < 2M (any stationary 
shell must necessarily be spacelike), as expected. Further details, including exact solutions 
of Eq. ( |3.17| ) , may be found in the original reference [0] . 
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IV. DUST SHELL WITH LEMAiTRE-TOLMAN-BONDI INTERIOR 



Let us now consider the case of a spherical dust shell matched to an interior LTB space- 



time [|T4|, and a Schwarzschild exterior. The interior metric reads 

ds^_ = -dt^ + F^(t, r)dr^ + X'^{t, r)dn^, (4.1) 



where F = X'/ a/1 — w{r), and AtcX"^ is the {t dependent) proper area of a spherical shell 
with coordinate radius r. For simplicity, we shall consider a class of LTB metrics given by a 
separable area radius function, X(t,r) = a(t)r. Included in this class are the LFRW cosmo- 
logical spacetimes, given by w{r) = Ar'^, where the constant A determines the geometry of 
the spatial section: hyperbolic, flat, or closed (spherical), for A = —1,0, or 1, respectively. 
The metric (O) is not in the standard spherical form ( |2.1| ), so the general formulae derived 



in Sec. II do not apply, but it is straightforward to compute Kab directly from its definition. 
In terms of the {x't} coordinates, the shell is given by 

t = T, r = Xc = const., (4.2) 

and the angular coordinates are, as before, trivially identified. The outward pointing unit 
normal to S is 

nt = ^^r. (4.3) 

a(r) 

where Wc = w{xc) = const., and continuity of the four-metric across S implies R{t) = 
a{T)xc. As before, we only need K^q-. 



K-=RVT^c, K+=R^1-^ + R\ (4.4) 

where M is the total gravitational mass appearing in the exterior Schwarzschild metric. The 
equation [K00] = —AnaR^ gives 

B^ = P-l + 2^1^ + §. (4.6) 

R R^ 

where /is = •AiraR^ and k = \/\ — Wc- This equation is formally identical to that for the 
two Schwarzschild shells [cf. Eq. (|3.9| )] provided we make the identifications k = k and 
/is = itl-e/2 = m+/(2/c). That is, a thin shell of proper mass ms moving in an interior 
Schwarzschild background with mass parameter m_ is equivalent to a thin shell of proper 
mass ms/2 moving in a LFRW background with total gravitational mass m_. The factor 



of 1/2 may be heuristically understood due to the presence of a continuous distribution of 
matter adjacent to the shell, instead of vacuum: half of the proper mass of the shell in the 
former case has already been accounted for in the inner side of the shell, in m_. 

Finally, we note that the usual boundary surface problem of matching a closed LFRW 
universe to a Schwarzschild exterior [|T^ is straightforwardly obtained from our construction 
by setting 

m-s = m+ = ^ M = m_, A = 1 ^ Wc = x^. 
The condition [Kab] = leads to a single non-trivial equation [cf. Eq. (^75|)] 



= 2^ - 1, (4.6) 

which is just the Friedmann equation for a closed universe {k = 1) with energy density 
/i = SM/iA-Kxl). 

V. SHELLS WITH PRESSURE 

We consider here the case of spherical shells with the following surface stress-energy tensor 
defined on them: 

Sab = (O- + p)UaUb + PJab, (5.1) 

where p is a tangential pressure, by virtue of spherical symmetry. Physically, this matter 
configuration can be realized by, e.g., 'counter- rotating' dust, where the individual dust 
particles are taken to follow geodesic orbits on the shell, with half of them orbiting along 
any given great circle with angular momentum per unit rest mass /, and the other half in 
the opposite direction with angular momentum per unit rest mass — /, such that the net 
angular momentum vanishes, thereby preserving spherical symmetry |16|. 



The junction conditions yield formally the same equation of motion as in the pressureless 



case [cf. Eq. (|3.8|) ], but, unlike in the latter case, the proper mass of the shell is no longer 
conserved. This is the crucial difference, and implies that the ratio m+ / ms is now an implicit 
function of R. The conservation equation (|2.14| ) gives 

cr + 2:^(^ + p) = 0. (5.2) 

Note that even the 'generalized' mass, ms = 47r(cr + 2p)R^ fails to be conserved in general. 
If the tangential pressure is constant, then is conserved, and the evolution equation 
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becomes identical to the one for the pressureless case, upon the substitution ijiy;, 
i.e., the constant pressure case is equivalent to the pressureless one, provided one includes 
pressure into the generalized mass definition (note, however, that even in this case ms 7^ 0, 
since there is always work done against pressure 'forces'). 

The conservation equation ( |5.2| ) can be rewritten as a first-order ODE for (j{R), by using 
a = {da/dR)R, as 

| + > + p)^0. (5.3) 

Upon introduction of the dimensionless function tp{R) = 27ia{R)R, the equation above can 
be rewritten as 

Now, one may either (i) specify an equation of state of the form p = p[a] and solve Eq. 
( p. 3D for cr(-R), or (ii) prescribe a functional form for ip{R), thereby fixing cr(i?), with p[a] 
following directly from Eq. ( ^.41 ). 

The equation of motion for the shell reads then 

i?2 + Y^ji^ = 0, (5.5) 

with the effective potential 



where 



a= {^] , b = M + m_. (5.7) 



The dynamics of the shell will obviously depend on ip and vice- versa, i.e., a given form of the 
potential — whereby the dynamics is uniquely determined — will constrain the allowed choices 
for ip, or equivalently, for the equation of state (although it will not uniquely determine 
them). The qualitative nature of the evolution is given by the shape of the effective potential: 
its zeros (if any), and first and second derivatives. From Eqs. ( |5.4| ) and (|5.6|) , we have then 



v 


= 2ip' 


V" 


= 2ij" 



. 2a b , . 



3 



37+0 -^IT,'^- (5-9) 



i?2^2 ' RJ R3 
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Now, in the presence of non-vanishing pressure it is possible to design the shape of V, so 
as to satisfy the necessary and sufficient conditions for a stable solution, R^: 



If the first inequality saturates, the solution is stationary. 

A. Non-existence of stable solutions with a strictly linear barotropic equation of 
state 

Proposition. There are no stable configurations with a strictly linear barotropic equation 
of state p = aa. 

Proof. Let us assume that p = acr, with a G (—1, 1), where the upper and lower hmits 
are imposed by causality. Equation (|5.4|) integrates to ijj = y4/i?^+^", where A > is an 
integration constant, and we have then 




(5.10) 




b cA^ 



(5.11) 



(5.12) 



V 



ti 



Aa{Aa-l)—R 



Aa-2 





(5.13) 




(5.14) 



(5.15) 



Now, at the local extrema, from Eqs. ( |5.13| ) and ( p.l4| ), after some manipulation we obtain 








(5.16) 



The solution _Ro = X^ 



-l/(l+4a) 




is stable if and only if Vq > 0, i.e., if and only if either of the 



following holds: 
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If (i) holds, then a G (—1/4,0), and we must take Xq~ in ( |5.15| ) with the further condition 
that the discriminant is greater than unity; however, this requires a > or a < —1/2, both 
of which contradict — 1/4 < a < 0. Now, let us examine possibility (ii): we have either 
a > or a G (-1, -1/4). If a > 0, the condition Xl + A^{\ + 2a)/(2aa) < is violated. If 
a G (—1, —1/2), the former condition is also violated, so we are left with a G [—1/2, —1/4). 
But this implies (since a is negative) that < 0, and thus we must take Xq' in ( ^.15] ) 
subject to either a > or a < —1/2, both of which are incompatible with a G [—1/2, —1/4). 
We conclude, therefore, that there are no solutions with V = and V" > 0. □ 

It immediately follows that all momentarily static solutions — which exist whenever 
V = — are unstable, i.e., there are no stationary solutions for shells with a strictly lin- 
ear barotropic equation of state. 



B. Energy conditions and stability criteria 



From Eq. ( |5.9|) , together with the requirements V = V = 0, after some algebra, the 
necessary and sufficient condition for a stable solution (V" > 0) may be written as 



a 



2 



+ 



> 



^'2 



+ 



Now, from Eqs. 



R^^P^ \^ r) R^ 2' 
we have 



(5.17) 



(5.18) 



where Cg = ^Jdp/da is the local sound speed, if one regards p as a hydrostatic pressure. 
The stability condition ( p.l7| ) reads then 

a 



{a+p)(l + 2Cl) ^ Q 
a + 2p ^ 



L27rV2/?6 



P(0 + 62 



where <^ = p/a, P(^) = 1 + 2<; + 61;^, and the 9j are strictly positive terms: 



62 



27rV2i^5 
a + 2-K%a^R^' 
^ + 27r'{a + 2pY. 



(5.19) 

(5.20) 
(5.21) 



The quadratic P(<j) has an absolute positive minimum at <j = —1/6, and therefore the 
right-hand-side of (|5.19| ) is strictly positive. 
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Now, in terms of the energy density and principal pressures, the energy conditions for 
this system are: 

Null energy condition (NEC): a + p > 

Dominant energy condition (DEC): cr > 0, cr > |p| 

Weak energy condition (WEC): a>0, a + p>0 

Strong energy condition (SEC): a + p>0, a + 2p>0 

The DEC imphes the WEC, which in turn imphes the NEC, and the SEC also implies the 
NEC; the energy conditions are otherwise independent. Comparison between these energy 
conditions and Eqs. ( |5.4| ) and shows that the latter provide a clear relation between 

the former and the matter content of the shell. Accordingly, one expects the stability, or 
lack thereof, of the shell to depend crucially on the energy conditions. 

For a strictly linear barotropic equation of state, condition (|5.19|) implies a + p > 0. 
Clearly, if the WEC is violated, there are no stable configurations. (In fact, as shown in 
the previous subsection, even if the WEC is satisfied, there are no stable solutions for such 
equations of state). 

For arbitrary equations of state, one may distinguish two cases, depending on whether 
is non-negative. Most macroscopic forms of matter obey > [including matter opaque to 
sound {C^ = 0), such as tenuous plasmas where electromagnetic radiation below the plasma 
cut-off frequency does not propagate], whereas certain fundamental fields, such as massive 
scalar fields, exhibit < 0. We shall assume throughout that cr > 0, remarking, however, 
that matter with negative energy density has been considered from a purely theoretical 
standpoint, e.g., in attempts to construct stable wormhole solutions ||T^. We skip trivial 
algebra and present the relation between energy conditions and stability in Table I below. 

VI. TRAPPED SURFACES 

Let S be any given compact spatial two-surface in the spacetime, and 6± be the ex- 
pansions in the future-pointing null directions orthogonal to E. The latter is said to be a 
trapped surface if 6_|_6_ > [0. Trapped surfaces signal thus the boundary of a region 
where any initially expanding null congruence begins to converge; clearly, they define regions 
of 'no escape'. The limiting case where expansion vanishes along future-oriented normal null 
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WEC SEC Unstable 



> 


Yes 


No 


Yes 


> 


No 


No 


Possibly 


> 


Yes 


Yes 


Possibly 


< 


Yes 


No 


-1/2 < C2 < 


< 


No 


No 


-1 < C2 < -1/2 


< 


Yes 


Yes 


-1 < C2 < -1/2 



TABLE I: Stability criteria and energy conditions. 

directions is referred to as an outer marginally trapped surface (OMTS). It is well known 
that the existence of trapped surfaces implies that of OMTS |jl9|. In spherical symmetry, 
the apparent horizon (AH) — the outer boundary of a compact trapped region — is an OMTS. 

Let us consider the case of an arbitrary shell S matched to two arbitrary spherical 
spacetimes, as described in Sec. II. Take the vacuum metric to be of the form (|2.1|) , with 
Qrr = '^1 F±{t, r), without imposing a priori any sign constraints on F. Since the four-metric 
g'^j^ is continuous across E and one needs derivatives of null vector fields orthogonal to E, it 
suffices to have g'^^^ to compute 9 along ingoing and outgoing null directions. Our surface 
S is that of a spherical shell with proper area radius r = R{t) (which provides a coordinate 
invariant definition, and, in addition, allows for measurements by an external observer, who 
can measure proper circumferences, but not coordinate radii). It is convenient to introduce 
null coordinates {u,v} defined by 

1 1 
du = hdt i=dr, dv = hdt H i=dr, (6.1) 



from which it follows that 

where the '±' sign denotes evaluation along outgoing and ingoing null directions, respectively. 
The condition for trapped surfaces reads then [where Eq. ( p.6|) was used] 

e+0_ = : F>0. (6.3) 

The AH is given by the limit O+0_ 0, which requires F — > 0. For the case of a boundary 
surface, [Kab] = 0, the time coordinate is globally defined, and thus F_ = everywhere, 
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FIG. 1: Apparent (AH) and event horizon (EH) evolution for a thin shell S with active gravitational 
mass m+ imploding in an interior Schwarzschild background with mass parameter m_. 

thereby uniquely defining a continuous curve tah('"), via F± = 0. In the presence of a 
thin shell, however, the functional forms of F_ and are different, and therefore the AH 
curve (implicitly defined by two different curves in A^±, by F± = 0, respectively) may not 
necessarily be continuous. 

To illustrate this, let us consider the case of a single shell with active gravitational mass 
m_|_ moving in an interior Schwarzschild background with mass parameter m_. This case 
is of sufficient generality, since, as shown in the preceding sections, the vacuum- immersed 
single shell well as that of a shell collapsing onto an interior LTB solution, can be 

trivially obtained from it. We have then 



1 



2m_ 



0, F^ 



0, 



(6.4) 



which defines the AH as 

tah(r) = 2m+Hi{r - R) + 2m_, (6.5) 

where Hi{x) is the unit Heaviside step- function {Hi = for x < 0, and Hi = 1, for a; > 0). 
The evolution of the AH is shown in Fig. |l]. 

The AH coincides initially with the EH, since the inner spacetime is just static 
Schwarzschild, but jumps out dis continuously at t = t*, when the thin matter shell S crosses 
the initial Schwarzschild radius. The EH is determined by the entire future development 
(along null generators) of the spacetime, and is therefore continuous, unlike the AH, which 
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is defined locally in terms of the gravitational mass interior to a certain radius, on a given 
spacelike hypersurface. 



VII. SHELL CROSSING 



In spherical gravitational collapse, one may distinguish between two different types of sin- 
gularities, given by the location and relative motion of neighboring shells. A shell focusing 
singularity is said to occur at the center of symmetry, when the proper area radius vanishes 
therein, leading to the blow-up of all curvature invariants; these singularities are gravita- 
tionally strong (any volume form defined along the Jacobi fields vanishes at the singularity), 
and the spacetime is geodesically incomplete pO|]. A shell crossing singularity is said to 
occur when neighboring shells cross each other at finite (comoving) radius, thereby leading 
to two-dimensional caustics, where the energy density diverges, and the metric becomes 



singular In LTB collapse, shell crossing occurs for spatially inhomogeneous density 

distributions, whenever the proper time for collapse of a given shell is a monotonically de- 
creasing function of the comoving radius. For the charged dust case, shell crossing has been 



shown to be inevitable, both for asymptotically flat p2|, and asymptotically de Sitter ||23 



spacetimes. Physically, shell crossings signal the intersection of matter flow lines at a given 
spacelike surface, to the future of which the model therefore breaks down. Such shell crossing 



singularities have been shown to be gravitationally weak p4[, and, in one particular LTB 
case, an analytical metric continuation was found in a distributional sense p5|. Recently, 
Szekeres and Lun ||2^ showed that it is always possible to find a coordinate transformation 
that renders the metric (but not C^) at the LTB shell crossing singularity. Their result 
adds considerably to our understanding of shell crossing in LTB collapse, but there remains 
to be shown that such transformation always exists (and is presumably unique) for a 
generic spherical metric. 



A. Conditions for shell crossing 

One might hope that a detailed analysis of the dynamics of individual thin shells could 
shed some light on the shell crossing process in continuous matter models. Thus motivated, 
we shall consider the case of two Schwarzschild shells, as described in Sec. III. Take the 
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shells to have initial radii R^, and dynamics given by 

+ V^{R^) = 0, (7.1) 

where 



2 



a_ = 2k\Jh- = rri-, 

a+ = 2{k^/hl + m_) = M + m_. 

For simplicity, and without loss of generality, we shall consider the marginally bound case, 
k = 1. In this case, the dynamical equation (|7.1| ) admits an analytical solution R^{t), which 
is explicitly given in Appendix A. The proper time for collapse is given by solving R{tc) = 0, 
which yields 



2 /„,,/o ,/o /463 3bR. 



?2 



where the '±' subscripts have been omitted for clarity. For given gravitational masses m± 
(whereby the corresponding a±, h± are determined), the proper time for collapse is a function 
of the initial radius R^. Now, shell crossing occurs iff, for Rq > Rq , we have 

rc(i?o+) < n{Ro)- (7.4) 

For given m±, it is straightforward to numerically solve for condition ( |7.4| ), thereby obtaining 
a region of the two-dimensional parameter space — given by the half-plane {Rq G (0, +cx3)} x 
{Rq > Rq} — in which ( |7.4| ) is satisfied. One can actually obtain an analytical criterion for 
shell crossing, by restricting our attention to neighboring shells, defined by 

i?+ = i?o(l + 0, (7.5) 
with < ^ ^ 1. For shells thus defined, we can Taylor expand 

re(i?+) = r,(i?o-) + ^_ e + 0{e). (7.6) 



"•0 



Condition (fT^ ) reads then 

^"^""'^ < 0. (7.7) 



dR, 



/ R- 
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FIG. 2: Gravitational collapse for two neighboring shells with data m+ = m_ = 1, Rq = 1/7, 
Rq = i^Q (1 + 10^^). The condition Rq < Rq < R^: (where i?* = 1/6) is obeyed, and hence the two 
shells cross before collapsing to zero radius. Note that both shells are initially inside the respective 
Schwarzschild radii. 

From Eq. ( [7.3|) , after a little algebra, this gives 



where R^, = 26_/a_ = m_/2. We note that this (necessary and sufficient) condition for 
shell crossing for neighboring shells requires the inner shell to be inside its Schwarzschild 
radius. From the definition of neighboring shells, it follows that the outer shell is also inside 
its Schwarschild radius: Take the limiting case Rq = m_/2, and set Rq = m_/2 + 6, 
with 5 ~ 0{mZ^)- Suppose that Rq is outside its Schwarzschild radius, i.e., -Rq > R^^^ = 
2(m_|. + m). This implies 



which cannot be, since m_|_ > and 6 ~ 0{mZ^); hence, Rq < R^^^. Note that the converse 
is not necessarily true: one may have -Rq > R^^^ and -Rq < -Rs^h) provided m+ < C(ml^). 

The case of non-neighboring shells [i.e., Rq — Rq ~ 0{Rq)] is qualitatively analogous: 
shell crossing only occurs if either (i) both shells are "sufficiently" inside their respective 
Schwarzschild radii, or (ii) the inner shell is initially untrapped, and the outer shell "suffi- 
ciently" trapped [0]. Illustrative examples are shown below. 



Ro e (0, R,) 



(7.8) 



4m+ + 3m_ < 26, 
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R 




FIG. 3: Gravitational collapse for two shells with data m+ = m_ = 1, /R^^-^ = /Rtch ~ 
Both shells are initially outside their Schwarzschild radii, and, as such, fail to obey the shell crossing. 




tau 



FIG. 4: Gravitational collapse for two shells with data m+ = m_ = 1, Ro/Rgch ~ 

Rq /Rgch ~ '^^^ outer shell is initially inside its Schwarzschild radius, but not the inner one. 

In this case, the shells cross at finite radius. 

B. Discrete shell model vs. LTB collapse 

In a multi-thin-shell model, the metric is at least at the shell crossing by construction, 
since [Qfiu] = for all times. This is in apparent contradiction with the well-known result for 



LTB collapse, where the comoving metric becomes singular {grr 0) at shell crossings [24 



Whilst such a metric singularity may be removed by a coordinate transformation via the 
prescription of Szekeres and Lun |26|, one is still left with the fact that scalar quantities 
such as the energy density remain divergent at LTB shell crossings (in fact, tidal forces 
are still infinite thereon), whereas they are finite in the discrete shell case [in particular, 
a = m/(47ri?g^) < oo]. Accordingly, the former cannot be taken as the continuous limit of 
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the latter. 

The inequivalence — as far as shell crossings are concerned — between these two models 
also arises from the following observation: In LTB collapse, individual shells follow geodesic 
motion in the four- dimensional spacetime due to the lack of any external forces (collapse is 
pressureless and uncharged), but this does not happen in the present case; within each shell, 
particle motion is geodesic, but, from the viewpoint of the four- dimensional spacetime, the 
shell is accelerated: 



u^'VuU^ = R 
h 



+ {R + 



R h' , 



K + 0, 



u 



b{3) 



(7.9) 
(7.10) 



We remark that the shell crossing process cannot be determined from the field equations 
alone, as a matter of principle. Once two (or more) shells cross, one must specify a priori 
the type of interaction that takes place (which depends on the microphysics of the model 
at hand), in order to follow the dynamics into the causal future of the shell cross surface. 
Only in special cases, where a purely gravitational gravitational interaction (whence the rest 
mass of each shell is conserved) between different shells is assumed, can one determine the 
evolution of the system beyond the shell cross. For null surfaces, we have the well-known 



DTR (Dray-t'Hooft-Redmount) relations which relate the mass and momenta of the 
different regions before and after shell crossing. Such relations were subsequently generalized 
to the timelike case (of massive spherical shells) by Nunez, Oliveira, and Salim p8 . 



VIII. CONCLUDING REMARKS 



We have analyzed in some detail the dynamics, stability, trapped surfaces, and shell 
crossing for spherical thin shells in vacuum. By focusing on the dynamics of the exterior 
shell in a vacuum two-shell model, the problem reduces to that of a single shell moving 
freely (albeit not geodesically) in a Schwarzschild background with different active mass 
parameters on each side. 

We have shown that there can be no stable solutions with a strictly linear barotropic 
equation of state, regardless of any energy conditions satisfied by such matter fields. Whereas 
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this rules out several idealized forms of matter as candidates for stable shell configurations, 
it still leaves ample room for stability with more general equations of state. For instance, by 
considering a quasilinear equation of state (which is a good approximation to any equation of 
state near equilibrium, if the latter exists), Brady, Louko, and Poison showed that stable 



static shells with an interior Schwarzschild background can exist, provided certain conditions 
are satisfied, which turn out to supersede energy conditions necessary for the existence 
(but not necessarily stability) of such shells |^. For matter configurations with > 



obeying the WEC, but otherwise arbitrary, we showed here that the SEC is a necessary 
and sufficient condition for the existence of stable shells around a classical Schwarzschild 
black hole, thereby producing a very simple and useful test for the stability of any matter 
content. For matter with G (—1/2,0) obeying the WEC, the SEC is also necessary 
and sufficient for stability. For G (—1, —1/2) the analysis cannot be carried out further 
without specifying an equation of state, to evaluate inequality ( ^.19|) . However, the results 



of for static shells with < suggest that there are no stable configurations, even 
dropping the staticity assumption. 

The collapse of a thin matter shell onto a Schwarzschild black hole introduces a step- 
function-type discontinuity in the apparent horizon curve, which occurs when the shell 
crosses the initial Schwarzschild radius, whence the jump equals the (active) gravitational 
mass of the shell. 

Neighboring dust shells were shown to cross whenever the inner shell is sufficiently inside 
its Schwarzschild radius, specifically, < R^^^il^- Shell crossings occur in a multi-dust- 
shell case, just as they do in the continuous LTB dust case, but individual shells move 
geodesically in the latter, whereas they are accelerated in the former. This, together with 
the fact that the energy density remains finite in the discrete case, but diverges in the 
LTB one, implies that the multi-shell case cannot be taken as the discrete analogue of LTB 
collapse, insofar as shell crossings are concerned. 
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APPENDIX A: MARGINALLY BOUND COLLAPSE OF TWO DUST SHELLS 



Formally, we need to solve the equation 



^+a/| + ;^ = 0, (Ai) 



with the initial condition R{0) = Rq. The positive real solution is 



where 



X (-1663 + 4a^E' - Ua'/'Er + OaV^) (A3) 
a"* a 



E = [4--3-^Rl + Rl] . (A4) 



When 6 = (i.e., m_ = 0), corresponding to a single shell m+ in vacuum, the above 
solution reduces to the well-known marginally bound case of LTB collapse: 

R(r) = (5) - yaK"r + ar') . (A5) 
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